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Discovery of a New Supernova Remnant in the Direction of 
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Abstract 

We discovered a middle-aged supernova remnant in the vicinity of G69. 7+1.0 using the ASCA satelHte. 
G69. 7+1.0 was identified in the 2.7 GHz survey and classified as a shell-type SNR with a diameter of 16'. 
During the ROSAT all-sky survey, the X-ray emission was detected in the direction of G69. 7+1.0. However, 
it extends beyond the radio shell, and an X-ray bright region was located outside of the radio shell. A 
spectral study with the ASCA and ROSAT shows a thin thermal plasma with an electron temperature 
of ~ 0.4 keV. There is no significant variation of the spectral parameters over the field of view, except 
for the lower column density of the eastern part. We also found a large shell structure which surrounds 
the X-ray bright region in both optical and radio images. We suggest that the observed X-ray emission 
is associated with the large optical and radio shell, and that they are part of a new SNR, different from 
the radio SNR G69. 7+1.0, which we have named AX J2001+3235 or G69.4+1.2. The large shell and the 
electron temperature of ~ 0.4 keV indicate that AX J2001+3235 is an evolved SNR. From a comparison 
with the column density of CTB 80 (G69.0+2.7), we estimate that the distance of the SNR is about 2.5 kpc. 
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1. Introduction 

Most galactic supernova remnants (SNRs) have been 
discovered and classified by observations in the radio 
band. Green's present catalog (1998 September version; 
Green 1998) lists 220 SNRs. These remnants are clas- 
sified into three types: "shell- type", "fiUed-center" , and 
"composite". Except for fiUed-center SNRs, they have 
radio shells. While some SNRs classified as shell- type in 
radio also have shells in X-rays, others show a centrally- 
brightened X-ray morphology. 

SNRs are also classified according to their age. While 
Cassiopeia A and Tycho's SNR, which are ~ 300 yr 
(Fabian et al. 1980) and 430 yr old (Reid et al. 1982), 
are classified as young shell-type remnants, the Cygnus 
Loop, Vela, and Puppis SNRs with ages of between 4000 
and 20000 yr (Ku et al. 1984; Gorenstein et al. 1974), 
are classified as evolved, middle-aged remnants. These 
evolved remnants typically appear as large radio shells 
containing many optical filaments. The radio emission 
is synchrotron radiation, which occurs when relativistic 
electrons move through a magnetic field, and the radio 
shell is considered to trace the blast wave. The X-ray 
emission comes from a thin thermal plasma heated either 
by the blast wave or by the reverse shock. Therefore, we 
can recognize evolved SNRs as X-ray emission well within 



a large radio shell (Seward et al. 1990; Green 1998). 

The SNR G69. 7+1.0 was discovered in the 2.7 GHz 
survey. The radio spectrum has an energy index, a {Si, ~ 
), of 0.8 (Reich et al. 1988) and it is classified as a 
shell- type SNR with a diameter of 16'. A S-D relation 
suggests the distance to G69. 7+1.0 to be 14.4 kpc (Case, 
Bhattacharya 1998). 

In X-rays, the ROSAT all-sky survey discovered a 
bright X-ray arc in the direction of G69.7+1.0 (Asaoka 
et al. 1996). It was observed in 1993 October with the 
ROSAT position-sensitive proportional counter (PSPC) 
in the pointing mode. The spectral study of this X-ray 
arc has not been made yet, which can clarify its nature. 
In this paper, we present the results obtained with the 
ASCA sateUite combined with the ROSAT PSPC data. 

2. Observations and Data Analysis 

Observations with the ASCA satellite (Tanaka et al. 
1994) in the direction of G69. 7+1.0 were performed on 
1997 November 20-21. We selected 3 positions, listed in 
table 1, in order to cover the entire emission based on the 
ROSAT observation. ASCA has four detectors mounted 
in the focal plane of its telescope. Two of the detectors 
are solid-state imaging spectrometers (SIS and SIS 1) 
(Yamashita et al. 1997), which were working in the 2- 
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Table 1. Summary of ASCA observations. 



Sequence 


Field center 




EflFective exposure 


SIS CCD chip* 


number 


a(J2000) 


J(J2000) 


[GIS/SIS] (ks) 




5502700 


20*'03"12^ 


32°43'56" 


7/10 


S0C1/S1C3 










S0C2/S1C0 


5502710 


20^00^51^ 


32°49'25" 


8/9 


S0C1/S1C3 










S0C2/S1CG 


5502720 


20''01'"34^ 


32°30'55" 


11/11 


S0C0/S1C2 










S0C1/S1C3 



* "S0C1/S1C3" is an abbreviation for "SIS0-chipl/SISl-cliip3". 



ROSAT PSPC 0.1-2.4 keV image 



+ 33=20' 



+ 33=00' 

o 
o 
o 

Z 

+ 32=40' 2 
H 
< 



'+:32S20''^5 



+,32=0.Q' 

20''04"' 20''02'" 2f}%0"" 

RIGHT ASCENSION 03MOt 




Fig. 1. ROSAT PSPC image witli tlic FOVs of tlic ASCA observation. The gray-scale bar is in units of counts s~^ arcmin"'^. 
The circles and the squares represent the CIS and SIS FOVs, respectively. The cross shows the position of G69. 7+1.0 
obtained with the radio observation (Green 1998). 



CCD Faint mode. The other two XRT are equipped with 
gas-imaging spectrometers (GIS 2 and GIS 3) (Ohashi et 
al. 1996; Makishima et al. 1996), which were working in 
the PH mode. While the GISs have a field of view (FOV) 
of 50' in diameter, the SISs have a rectangular 11' x 22' 
FOV in the 2-CCD mode. Figure 1 shows the ASCA 
FOV superposed on the ROSAT PSPC image. Detailed 
information concerning the observation are also listed in 
table 1. 



We rejected the data at the South Atlantic Anomaly 

(SAA), the high-background regions with the geomag- 
netic cut-off rigidities of < 6 GeV/c. We also excluded 
the data at low-elevation angles from the Earth rim. The 
data screening criteria are < 5° from the night earth rim, 
< 25° from the bright earth rim for the GIS, < 45° for 
the SIS 0, < 20° for the SIS 1. As for the GIS data, 
we applied a "flare-cut" to maximize the signal-to-noise 
ratio described in Y. Ishisaki et al. (1997 ASCA News 5, 
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Fig. 2. GIS mosaic images observed in the 0.7-1.5 kcV band(a) and in tlie 1.5-5 keV band(b). The scale bars are in units 
of counts s~^ arcmin"^. The rectangles and circles in figure 2a indicate the regions from which the spectra displayed in 
figure 3 were extracted. 



26). SIS data were acquired in the 2-CCD mode in 1997, 
and hence the residual dark distribution (RDD) wotdd 
become a problem. Therefore, we applied the RDD cor- 
rection to the SIS data. Then, they were corrected for 
any dark-frame error, the echo effect and any charge- 
transfer inefficiency (T. Dotani et al. 1995 ASCA News 
3, 25, 1997 ASCA News 5, 14). Hot and flickering pixels 
were also removed. 

The background consists of the non- X-ray background 
and the cosmic X-ray background. For the GIS, the non- 
X-ray background was reproduced with a method intro- 
duced by Ishisaki (1996). The cosmic X-ray background 
data were extracted from the Large Sky Survey (Ueda 
et al. 1999) data observed during the ASCA PV phase. 
We subtracted these background data from the GIS data. 
For the SIS, we subtracted the blank-sky data as back- 
ground. We estimated the contribution of the galac- 
tic ridge emission to be negligible (Yamauchi, Koyama 
1993). In a spectral analysis, we binned the data into 
groups having a minimum of 20 counts per bin, and sub- 
tracted the background spectrum, which was grouped to 
match the source data. 

As for the ROSAT PSPC observation, we retrieved 
the ROSAT archival data through the HEASARC/GSFC 
Online Service. The sequence ID of the data used here 
is RP500323. The pointing position was (a, (5)(J2000) = 
(20'"02"^40?8, 32°43'48."0) and the net exposure time was 
12 ks. The spectra were rebinned as well as those of the 
ASCA data. The background spectrum was determined 



from a circular region of radius 34' at the center of the 
PSPC FOV, excluding both the X-ray arc and the point 
sources. 

3. Results 

3.1. X-Ray Image 

We generated 
exposure-corrected, background-subtracted images of the 
GIS. Figure 2 shows GIS mosaic images both in the soft 
band (0.7-1.5 keV) and in the hard band (1.5-5 keV). 
The soft band image (figure 2a) shows an extended X- 
ray emission which is similar to the extended source "A" 
seen in the ROSAT PSPC image (figure 1). The other 
brightest sources in figure 1 are point sources "B" and 
"C" . Source C is missing in the ASCA image. Since the 
PSPC count rate of source C is ~ 30% of that of source 
B, the expected fluxes of source C for the GIS and SIS 
are considered to be below the threshold. The hard band 
image (figure 2b) shows only two sources: the northern 
source is a point source from the ASCA PSF of the XRT, 
whereas the southern source is surely an extended source. 

The northern point source in figures 2a and b was 
identical to source B seen in the ROSAT PSPC ob- 
servation. Source B corresponds to the ROSAT All- 
Sky Survey Bright Source Catalogue (RASS-BSC) source 
IRXS J200047.0-h325708 at position (a, (5)(J2000) = 
(20'^00'"47f0, 32°57'08."5), and it is within 20" from a 
K2 star, HD 189806 (TYC 2674-1400-1), to which the 
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Table 2. Best-fit parameters of the thin thermal model (VMEKAL model).* 



Parameter 


Bright part 


Eastern part 


Northern part 


Radio shell 


Entire^ 


TVh (10^1 cm-2) 
fcTc (keV) 


1.9±0.4 
0.45±0.03 


0.7±0.1 

n 46+" °^ 

^■^"^-0.02 


■^•3-0.6 

40+° °'' 

^•^^-0.03 


^■^-0.3 

n 47+0-08 

^•^'-0.13 


-■^■"-0.4 


Abundance (relative to 


cosmic value) 










Si 


0.3±0.2 


^•^-0.3 


0.8±0.4 


< 2.0 


0.6±0.2 


Fe 


0.07±0.03 


0.07±0.03 


< 0.09 


16+°°^ 

^■-'^"-0.10 


0.13±0.04 


xVd.o.f 


311.8/269 


315.7/274 


288.1/261 


98.2/132 


685.6/445 



* Errors on parameters are computed using A^^ = 2.7. 
I Dashed elliptic region shown in figure 2a. 



distance is 56 pc (Voges et al. 1999). 

We extracted the spectra of this source from a 3' radius 
circle for the GISs, SISs and ROSAT PSPC. In order to 
subtract the component of extended emission from the 
data, the background spectra were determined from a 
4'-8' radius annular region for the GISs and PSPC, and 
from the FOV of the CCD chip S0C1/S1C3 excluding a 
4' radius circle of the source for the SISs. Since it has 
been reported that SIS data tend to infer a higher column 
density (ASCA Guest Observer Facility's home page at 
http://heasarc.gsfc.nasa.gov/docs/asca/caLprobs.htm ), 
we used SIS data above 0.8 keV. 

We performed a combined-fit using five data sets, 
GIS 2, GIS 3, SIS 0, SIS 1, and PSPC data. Since we 
could expect coronal X-ray emission from the K2 star, we 
fitted them with a thin thermal model (MEKAL model, 
Mewe et al. 1985; Liedahl et al. 1995) of cosmic abun- 
dances (Anders, Grevesse 1989). We used the model of 
Morrison and McCammon (1983) for the interstellar ab- 
sorption. The parameter values that provide a minimum 
were taken as "best-fit" parameters, and errors on 
parameters were computed using Ax^ = 2.7. Although 
the MEKAL model gave only a poor fit with x^/d.o.f 
= 198.8/99, we obtained the best-fit values for the ab- 
sorption feature, TVh < 3 x 10^* cm~^ and the electron 
temperature, kTc = 0.28±0.01 keV, respectively. The ob- 
tained value for kTe is typical for the temperature of the 
coronal plasma of late-type stars (Schmitt et al. 1990). 

3.2. X-Ray Spectra 

We extracted the spectra from those regions where the 
SIS observed. There are three regions: an eastern part, 
a northern part and a bright part. As for the data of 
the eastern part and the northern part, we excluded the 
regions (marked with circles inside the SIS rectangles in 
figure 2a) within 4' of sources B and C seen in the ROSAT 
PSPC image. In addition to these three regions, we ob- 
tained the spectrum of the region which corresponds to 
the bright region in radio, shown in figure 2a as a cir- 



cle. The circle has a diameter of 10'; since it is out of 
the SIS FOVs, we used the GIS data for this spectrum. 
The ASCA and ROSAT PSPC spectra of aU regions are 
shown in figure 3. 

We tried to fit spectra with a MEKAL model with vari- 
able metallicity, modified by the interstellar absorption. 
In this case, there are residuals at around 0.8-1.3 keV 
and 1.8 keV. Therefore, we used a VMEKAL model and 
allowed the abundances of Si and Fe to vary during the 
fitting. The abundances of other elements were set to 
be cosmic values. The best-fit curves are shown by the 
solid lines in figure 3, and table 2 summarizes the best- 
fit parameters. Although the eastern part has a slightly 
lower A'h value than other regions, there is no significant 
variation in the electron temperature over our FOV . 

In figure 3, the remaining residuals above 2 keV can 
be seen at the spectrum of the bright part, which may 
indicate an additional non-thermal component. We then 
fitted the spectrum with a two-component model con- 
sisting of a thin thermal model and a power-law model. 
Since we could not determine all of the parameters when 
the electron temperature and the abundances of Si and 
Fe are free parameters, we fixed these parameters to the 
values which we obtained with the fitting using only a 
thin thermal model. However, this attempt could not 
produce a significant improvement to the fit. 



3.3. Timing Analysis 

We searched for pulsations using GIS data within 4' 
from the peak of the bright region. We used both high 
and medium bit-rate data. After barycentering the pho- 
ton times of arrival, the power spectral density was cal- 
culated by the fast Fourier-transformation method. The 
Nyquist frequency was set at 8 Hz. We found no signif- 
icant pulsation. Other regions in our FOV do not have 
sufficient counting statistics to perform a reliable timing 
analysis. 
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Fig. 3. Top panels show the X-ray spectra obtained with the ASCA and ROSAT PSPC in each case. Best-fit curves are shown 
with solid lines, (a)-(c) From top to bottom: PSPC (xlOO), CIS 3 (xlO), CIS 2, SIS 1 (xO.l), SIS (xO.Ol), (d) PSPC 
(xlOO), CIS 3 (xlO), CIS 2. The lower panels show the residuals of the fits for only CIS 2 and SIS besides the PSPC. 



4. Discussion 

Wc compared the X-ray emission with the radio emis- 
sion in the direction of G69. 7-1-1.0. Figure 4 shows the 
ROSAT PSPC image superimposed on the 4850 MHz 
radio contours. The radio map was obtained through 
the Skyview supported by HEASARC/GSFC. The X-ray 
emission observed by the PSPC extends beyond the radio 
shell of G69.7-hl.O, which is centered at {a, (5)(J2000) = 
(20'^02'"42^ 32°43'28"). The X-ray peak occurs roughly 
at (a, 5)(J2000) = (20'^01™10% 32°35'), well outside the 
radio shell, at 25' from its center. Usually, the radio shell 
of an SNR appears to correspond to, or surround, its 
X-ray bright region (Seward 1990; Green 1998). There- 



fore, it is unlikely that the bright part of the X-ray emis- 
sion is associated with G69. 7+1.0, but it is still prob- 
ably an SNR, because of its thermal emission and ex- 
tended morphology. We suggest that the observed X-ray 
emission may be a new SNR, different from the radio 
SNR, G69.7-hl.O, and have named it AX J2001-h3235, or 
G69.4-I-1.2. 

The radio map shows that a broken shell structure 

seems to surround the extended X-ray emission. This 
shell has a diameter of about 1°, as shown in figure 4. We 
also found a corresponding optical shell in the Digitized 
Sky Survey image (retrieved from HEASARC/GSFC, fig- 
ure 5). Even though additional observations would be 
required to establish whether the optical and radio emis- 
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Table 3. Parameters for evolved SNRs. 



SNR Distance kT^ Radius Et Age 

(kpc) (keV) (pc) (10^9 (10^ yr) 

AX J2001+3235 2^5 0.42^^;^^ 22 2.6+^ j 3.7 ±0.3 

lAA 0.42+g:g| 128 (2.1ib0.2) X 10^ 4.9 ±0.4 

CygnusLoop OJT 0^26^^ 5 T 

Vela SNR 0.25^ 0.11, 0.45*'" 17* 2 1* 

Puppis A 2.2§ 0.60§ 8.5§ 4 0.4§ 



* Ku et al. (1984) 

t Cha et al. (1999) 

X Bocchino et al. (1999) 

§ Gorcnstcin ct al. (1974) 

tt Bocchino et al. (1999) used a two-temperature emission model. 



sion is indeed an SNR, we note that this structure is 
remarkably similar to that observed in known SNR, as 
for instance, Cygnus Loop (Ku et al. 1984). Therefore, 
it cannot be excluded that this structure represents a 
shell of an evolved SNR. The X-ray emission with the 
temperature of about 0.4 keV supports that the SNR 
is in middle-aged. No association of the X-ray emission 
with this large shell has been reported so far. Asaoka 
et al. (1996) suggest that the X-ray emission is an SNR, 
and is a part of another shell with a radius of about 
3°, which is different from our conclusion. Our finding 
is that the 1° diameter shell at optical and radio wave- 
lengths may be related to the extended X-ray emission 
of AX J2001+3235. 

We note that the X-ray emission is elongated to the 
south-west direction where the radio and optical shells 
are broken. Both the elongation in X-rays and the break- 
out at optical and radio wavelengths may be due to 
the low density of the interstellar medium in the west- 
ern part. This fact confirms that the X-ray emission in 
the direction of G69. 7+1.0 is associated with an evolved 
SNR having both radio and optical shells, rather than 
G69.7±1.0. Follow-up radio and optical observations are 
needed to study more in detail the relationship between 
the X-ray emission and the 1° diameter shell. 

We have estimated the distance of AX J2001+3235 us- 
ing the measured column density. Wc extracted the CIS 
and PSPC spectra from the region shown as a dashed 
ellipse in figure 2a and fitted them with the same model 
as that described in subsection 3.2. The fitting result 
is shown in table 2. Although the poor fit suggests that 
the model is too simple, we employed this result to obtain 
the average SNR properties as a first-order estimate. Wc 
compared the column density of AX J2001±3235 with 
that of CTB 80 (G69.0+2.7) which is the SNR located 1° 



from AX J2001+3235. The distance to CTB 80 has been 

measured by various methods (dispersion measure of the 
central pulsar, kinematic estimates from the detected in- 
frared shell, and the optical reddening estimates), and 
2.5 kpc is commonly used (Blair et al. 1984; Case, Bhat- 
tacharya 1998). On the other hand, the best-fit value 
of the column density A^h of CTB 80 obtained with the 
ROSAT observation was 3.0 x lO^i cm'^ (Safi-Harb et 
al. 1995). The similar value of the column density of AX 
J2001+3235 A^H ~ 2 X 10^^ cm'^ suggests the similar 
distance of ^ 2.5 kpc. 

Employing the distance of 2.5 kpc, we calculated the 
parameters of AX J2001+3235. The fitting result of the 
spectra of a dashed elliptic region gives a lower limit 
of the thermal energy contained within the entire SNR. 
Prom the spectral fitting, we obtained kTg of 0.42 keV 
and EM/(47rrf2) ^ 4.0x10"!*' cm-5 where EM is the 
emission measure, / n^riB. dV , rie and riH are the mean 
number density of electron and hydrogen, and d is the 
distance to the SNR. Assuming that the object is an el- 
lipsoid with a major axis of OSl (~ 22 pc) and a minor 
axis of 024 (~ 10 pc), and the number density of hydro- 
gen is equal to that of electron, the derived mean number 
density of hydrogen nn is ~ 0.06 cm~'^. The thermal en- 
ergy Et is calculated to be 

Et ~ nnVkT^ 

where V is the volume of the calculated region. 

We compared the parameters of AX J2001±3235 with 
those of other evolved SNRs, the Cygnus Loop (Ku et 
al. 1984), Vela and Puppis SNRs (Gorenstein et al. 1974; 
Bocchino et al. 1999; Cha et al. 1999). The results are 
summarized in table 3. The age of AX J2001+3235 was 
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Fig. 4. ROSAT PSPC image of the vicinity of G69. 7+1.0 with an overlay of 4850 MHz survey contour map. The contour 
levels are 5, 10, 20, 40, 80, 160, 320 and 640 in units of mjy/beam. 



derived using the Sedov model with an explosion energy 
of 10'''^ erg. Wc can sec that all of the parameters (elec- 
tron temperature, radius, thermal energy, age) of AX 
J2001+3235 have typical values of evolved SNRs. We 
also show the parameters with the assumption of the 
distance of 14.4 kpc, which was obtained with the S- 
D relation for G69. 7+1.0. In this case, the size and the 
thermal energy of AX J2001+3235 arc much larger than 
other evolved SNRs. Therefore, a distance as large as 
14.4 kpc is very unlikely. 

5. Summary 

We have studied the X-ray emission of a large region 
around the known radio SNR G69.7+1.0, and found that 
the large extended X-ray source, AX J2001+3235, is not 
associated with G69. 7+1.0. This source has an elongated 
shape with a size of 1°, which is much bigger than that 
of G69.7+1.0. It is surrotmdcd by another large shell at 
both optical and radio wavelengths. The column density 
and the electron temperature obtained by the ASCA and 



ROSAT PSPC data are A^h ~ 2 x lO^^ cm'^ and kT^ 
- 0.4 kcV. The Ar value of AX J2001+3235 indicates 
that it is closer than 14.4 kpc, which is the distance to 
G69.7+1.0 derived from the S-D relation. We estimate 
the distance to AX J2001+3235 to be - 2.5 kpc. The 
radio and optical shells and the derived parameters of AX 
J2001+3235 using the fitting results of the X-ray spectra 
are consistent with the idea that AX J2001+3235 is an 
evolved, previously unknown, SNR. 

The authors are grateful to all the members of the 
ASCA team. We are also grateful to the referee, F. 
Bocchino, for his useful comments. K.Y. is supported 
by JSPS Research Fellowship for Young Scientists. 
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